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ABSTRACT 


Two significant contributions have been made during this research period in the 
research "Robust Control of Systems with Real Parameter Uncertainty and Unmodelled 
Dynamics" under NASA Research Grant NAG-1-1 102. They are: (1) a fast algorithm for 
computing the optimal H“ norm for the four-block, the two block, or the one-block optimal 
H~ optimization problem, and (2) a construction of an optimal H~ controller without 
numerical difficulty. 

In using GD (Glover and Doyle) or DGKF (Doyle, Glover, Khargonekar, and 
Francis) approach to solve the standard H“ optimization problem, the major computation 
burden is on the computation of the optimal H°° norm which required bisection search. In 
this research period, we developed a very fast iterative algorithm for this computation. Our 
algorithm was developed based on hyperbolic interpolations which is much faster than any 
existing algorithm. The lower bound of the parameter, y, in the H“ Riccati equation for 
solution existence is shown to be the be square root of the supremum over all frequencies 
of the maximum eigenvalue of a given transfer matrix which can be computed easily. The 
lower bound of y such that the H~ Riccati equation has positive semidefinite solution can be 
also obtained by hyperbolic interpolation search. Our another significant result in this 
research period is the elimination of the numerical difficulties arising in the construction of 
an optimal H“ controller by directly applying the Glover and Doyle's state-space formulas. 

With the fast iterative algorithm for the computation of the optimal H°° norm and the 
reliable construction of an optimal H“ controller, we are ready to apply these tools in the 
design of robust controllers for the systems with unmodelled uncertainties. These tools will 
be also very useful when we consider systems with structured uncertainties. 
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PROGRESS REPORT 


ROBUST CONTROL OF SYSTEMS WITH REAL PARAMETER 
UNCERTAINTY AND UNMODELLED DYNAMICS 

1. INTRODUCTION 

This document is the first-period progress report on the NASA supported research, 
"Robust Control of Systems with Real Parameter Uncertainty and Unmodelled Dynamics", 
(No. NAG- 1-1 102). The objective of the proposed research is to develop reliable and 
efficient algorithms for the computational problems arising in the design of robust optimal 
controllers for the systems with structured and unmodelled uncertainties and apply them to 
practical aerospace control systems. We are happy to report that we have obtained some 
significant results [1,2,3] in the solution of the four-block, the two block, or the one-block 
optimal H~ optimization problem. Meanwhile, we also have some primitive results [4,5] in 
the mixed fP/H” optimization problems which may have great impact on our future 
research. Furthermore, the work to be performed in the succeeding budget period has also 
been set on the right track. 

In the paper, "Iterative Computation of the Optimal H°° Norm By Two-Riccati- 
Equation Method" [1], an iterative algorithm for computing the optimal H“ norm by using 
DGKF two-Riccati-equation approach [6] was proposed. The celebrating two-Riccati- 
equation solution to a standard H°° control problem can be used to characterize all possible 
stabilizing optimal or suboptimal controllers if the optimal or suboptimal H 0 * norm is 
given. No efficient algorithm for computing the optimal or suboptimal H” norm was 
available in the literature before this paper was written. 

A numerical difficulty caused by the inversion of a singular matrix usually arises in 
using the DGKF state-space formulas to construct an optimal H“ controller. In the paper, 
"Design of an H~ Optimal Controller by DGKF State-Space Formulas" [2], we explained 
that the numerical difficulty originated from restricting the controller to be strictly proper 
and showed that the numerical difficulty can be easily eliminated if a proper controller with 
direct feedthrough is allowed. 

In the paper, Computation of the Optimal H°° Norm by H°° Riccati Equations and 
Hyperbolic Interpolations" [3], we proposed an extremely fast iterative algorithm for 
optimal H“ norm computation for a more general H“ optimization problem considered by 
Glover and Doyle [7]. The algorithm was developed based on hyperbolic interpolations 
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which is faster than any existing methods in the literature. The numerical difficulties arising 
in the construction of an optimal H“ controller by direct applying the Glover and Doyle's 
state-space formulas were also addressed in the paper. 

The papers, "Necessary and Sufficient Conditions for Mixed H 2 and H“ Optimal 
Control" [4] and "Design of a Suboptimal H°° Controller with H 2 and Bandwidth 
Constraints" [5], summarized our primitive results in the mixed H^/H” optimal control. 
These results prompt a new research problem which will be briefly described in section 3. 

In section 2 of this report, we will show the overall progress in this research 
period. In section 3, the work for future research will be briefly described. 


2. OVERALL PROGRESS 
2.1 Introduction 

Consider the system 


z(s) 

_y(s)_ 


G n (s) 

G 21 (s) 


g 12 (s) 

G 22 (s) 


v(s) 

u(s). 


G(s) 


v(s) 

|_u(s) 


(2.1- la) 


u(s) = K(s) y(s) (2.1-lb) 

where G n (s) 6 R(s) PlXmi , Gi 2 (s) e E(s) PlXm2 , G 2 i(s) e lR(s) P2Xmi , and G 22 (s) e 
R (s) P2Xm 2 R ( s fl»l 

is the set of pjxmj proper rational matrices with real coefficients. 
Recall that the standard H°° optimization problem is the problem of finding a proper 
controller K(s) such that the closed-loop system is internally stable and IID r i (G,K)ll eo is 

minimized where 

CT x (G,K) := G n + G 12 K (I - G^K)- 1 G 21 (2. 1 -2) 

is the transfer function of the closed-loop system from v to z and 

IIOII^ := sup a [O(jco)] (2.1-3) 

CO 

and °(‘) denotes the maximum singular value. 

Let 
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G(s) = 


(2.1-4) 


' A 

B, 

B„ " 


l 

2 

C l 

D 11 

D 12 

5 

D 21 

D 22_ 


be a minimal realization of G(s) and A e IR™ Here 



{ A, B, C, D } 


(2.1-5) 


implies a state-space realization and G(s) = D + C (si - A)' 1 B. 

According to the dimensions of G n (s), G 12 (s), G 21 (s), and G 22 (s), the standard 
H°° optimization problem has the following situations to consider: (a) p t > m 2 , p 2 < (b) 

Pi ^ m 2 , p 2 < nq; (c) Pi > m 2 , p 2 > nq; and (d) pj ^ m 2 , p 2 £ mj. Situation (a) is referred 
as the four-block H°° optimization problem. Situations (b) and (c) are referred as the two- 
block H°° optimization problem. Situation (d) is referred as the one-block H°° optimization 
problem. 

Recently, Glover and Doyle [7] presented a celebrating Riccati-equation type 
solution to the general H°° optimization problem (including four-block, two-block, and one- 
block problems). Glover and Doyle's approach characterizes all possible stabilizing 
suboptimal H°° controllers which order is not higher than that of the plant. In utilizing these 
Glover and Doyle's formulas to design an optimal (or suboptimal) H~ controller, the most 
computationally demanding work is the computation of the optimal (or suboptimal) H°° 
norm which requires iteration. Up to now in the literature, there is no efficient and 
systematic search algorithm available for this computation. 

In this report we will present a very fast iterative hyperbolic search algorithm for the 
computation of the optimal (or suboptimal) H“ norm. The optimum can occur in three 
cases. In case (1), the optimum occurs at the smallest y such that the two H~ Riccati 
equations have stabilizing solutions X and Y, i.e„ these X and Y happens to be positive 
semidefinite and p(XY) < y 2 . Case (2) occurs when Y=0 (or X=0) for all y and the optimal 
H~ norm is the smallest y such that X (or Y) is positive semidefinite. The most likely one to 
happen most of the time is case (3) in which the optimal H“ norm is the y such that the two 
H“ Riccati equations have positive semidefinite stabilizing solutions X and Y and p(XY) = 
y 2 where p(XY) is the spectral radius of XY. In subsection 2.3, the optimal H°° norm in 
case (1) is shown to be square root of the supremum over all frequencies of the maximum 
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eigenvalue of a given transfer matrix which can be computed easily [8,9]. The optimum in 
either case (2) or case (3) can be obtained by hyperbolic interpolation search which is 
presented in subsection 2.4. 

Except case (1) which does not happen very often, the Glover and Doyle's 
controller formulas [7] can not be directly used to construct an optimal controller since an 
inversion of a singular matrix will cause numerical difficulty. As mentioned in [6] and [10], 
a descriptor (or generalized state-space representation) version of the controller formulas 
can avoid this numerical difficulty. With slight rearrangement, the Glover and Doyle's 
controller formulas can be rewritten in a descriptor form. The descriptor form of an optimal 
controller can then be reduced to an at least one order less state-space representation which 
is addressed in subsection 2.5. 

For convenience, the Glover and Doyle's state-space formulas for suboptimal H~ 
solutions are listed in Subsection 2.2. In Subsections 2.3 and 2.4, algorithms to compute 
the optimal H~ norm in three cases are presented. The numerical difficulty in the 
construction of an optimal H“ controller is addressed in Subsection 2.5. Some illustrative 
examples are included in Subsection 2.6. 


2.2 Glover and Doyle's State-Space Formulas 

In [7], Glover and Doyle assume the realization of G(s) is given by (2.1-4) with the 
following assumptions. 

(i) (A, B 2 ) is stabilizable and (C^ A) is detectable. 

(ii) rank D 12 = m 2 , rank D 21 = p 2 . 


(iii) D i2 = 


> D 21 = [0 I], and D n is partitioned as 


^1111 ^1112 


L D 1121 D 1122 J 


with D 1122 e ]R P2Xm2 ) 
(iv) D 22 = 0 (this can be removed [7] ). 


jcol - A B 0 


(v) rank 


D 


= n + m^ V cue R. 


12 
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(vi) rank 


jcol - A 




= n+p 2 


V coelR. 


Define two Hamiltonian matrices as follows, • 


H„(y) := 


’ A 

0 


B 


T 

T 

_ 

T 

r '[dT. c , b t ] 

-C c 
L 11 

-a t 


L' C i D i.J 


(2.2- la) 


and 


UY) := 


A T 0 
-BjB’J' -a 


T 

- b , d .i 


R U 


D ., B I C ] 


(2.2- lb) 


where 


D,.-[ D 11 D 12 ] . D., = 


D 


n 


D 


21 


(2.2- lc) 


and 


R = D D, 

l* l* 


Vi mI 

1 

O 

Q 

II 

lOi 

1 

*— * 

0 

. 0 

0 

1 

• 1 *1 

. 0 OJ 


(2.2- Id) 


Then the following theorem shows an easy way to construct a suboptimal stabilizing 
controller such that liy ^(G,K)II^ < y where ^(G.K) is the closed-loop transfer matrix 

from v to z. 

Theorem 2.2-1: [7] 

There exists a stabilizing controller such that IIJT ^(G.K)!^ < y if and only if the 
following three conditions hold. 

(i) y > max ( a Pini ^ 1112 ]’ a ^1111 ^ 1121 ]) (2.2-2a) 

(ii) H^(y) e dom(Ric) and X(y) := Ric[H 00 (y)] > 0. 

J«(Y) 6 dom(Ric) and Y(y) := Ric[J 00 (Y)] > 0. (2.2-2b) 

(iii) P[X(y)Y(y)] <Y 2 - (2.2-2c) 

Moreover, when these conditions hold, one such controller is 
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K sub < s > = 


(2.2-3a) 


1 

>> 


e 

£> 


where 

^ = '^ 1121 ^ 1111 ^ ' D iii2"^ii22 • (2.2-3b) 

£= {F 2 -6(C 2 +F 12 )}Z (2.2-3c) 

6 = -H 2 + (B 2 +H 12 )6 (2.2-3d) 

A = A+HC+(B 2 +H 12 )£ (2.2-3e) 

Z := (I-y^YX)' 1 (2.2-3f) 

pT = [^n F^] = - (XB + c[ Dl .) R' 1 (2.2-3g) 

H = [ H i i H 12 H 2 ] = -(YC T + B 1 Dy i )R 1 (2.2-3h) 


and F n € R < m l-P2)™ p^ e ]rP 2“ p 2 € R m 2™ Hn e R nx( Pl -m 2 ^ Hw € E "xra 2j 
H 2 e ]R nxp2 . 

In the above theorem, condition (ii) means that there exist positive semi-definite 
solutions X and Y to the algebraic Riccati equations corresponding to the Hamiltonians 
FUy) and J«(y) respectively. Condition (iii) means that the spectral radius of XY is less 
than y 2 . 

The above theorem shows an easy state-space approach to construct a stabilizing 
suboptimal controller such that HIT ^(G.Kjll^ < y. The order of the suboptimal controller 

can be the same as that of the plant G(s). The major computation involved is the solution of 
two H“ Riccati equations which are easy to solve if solutions exist. 

The Glover and Doyle’s approach is a great break-through in the solution of H°° 
optimization problem. Theorem 2.2-1 can also be used to compute the optimal H“ norm 
and to construct an optimal H°° controller. Algorithms for computing the optimal H°° norm 
will be presented in subsection 2.4 and the construction of an optimal H~ controller will be 
addressed in subsection 2.5. Subsection 2.3 gives the lower bounds of y for the H°° Riccati 
equations to have solutions and to have positive semi-definite solutions respectively. 
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2.3 Lower Bounds for Solution Existence of H°° Riccati Equations 

As mentioned in the previous subsection, Theorem 2.2-1 can also be employed to 
compute the optimal H~ norm. The optimal H~ norm is the smallest y such that the three 

conditions in Theorem 2.2-1 are satisfied. The computation of the optimal H“ norm is 
closely related to the algebraic Riccati equations associated to the Hamiltonian matrices 
HJyjandJJy). 

For the algebraic Riccati equation associated to the Hamiltonian matrix H^fy) to 
have a unique stabilizing solution, H„,(y) must have stability and complementarity 
properties [6]. Stability property means that H^fy) has no eigenvalues on the jco-axis. With 
stability property, H^Xy) has n eigenvalues in Re s < 0 and n in Re s > 0. Let X.CHJ be an 
n-dimensional spectral space corresponding to eigenvalues in Re s < 0. Partitioning the 
matrix constructed from the basis vectors of X.OTJ, we have 



where X lt X 2 e R nxn . Complementarity property means that X t is nonsingular, or 
equivalently, the two subspaces 


X_(HJ, Im ° (2 - 3 ‘ 2) 

are complementary. When the Hamiltonian matrix H»(y) has stability and complementarity 

properties, X=X 2 X J 1 is the unique stabilizing solution of the algebraic Riccati equation 
associated to H^fy). 

Let a x be the smallest y such that the Hamiltonian matrix H M (y) has no jco-axis 
eigenvalues. Then for every y > a x except some isolated y at which X ( is singular, the 
algebraic Riccati equation associated to H»(y) will have a unique stabilizing solution. The 
smallest y such that the Hamiltonian H^fy) has no jco-axis eigenvalues is given by the 
following theorem. 

Theorem 2.3-1: 

The smallest y such that the Hamiltonian H^fy) has no jco-axis eigenvalues, denoted by 
a x , is the square root of 

S“P I a; , [I - o , 2 (0; 2 G 12 )‘ 1 g; 2 ]Q , , (jo) } (2.3-3) 


1 0 



where G n (jco) (G 12 (jco) resp.) is the conjugate transpose of G n (jco) (G 12 (jco) resp.) 
and Xmax means the maximum eigenvalue. 


Proof: 

Let 


r(s) = 


Gy i (-s)G n (s)-7 2 I ml 

Gj 2 (-s)G n (s) 


G T n (-s)G n (s) 


It is straightforward to show that 


(2.3-4) 


r ! (s) 


a-br’d^.c, 

-1 T 

-BR B 

-BR' 1 

- C I C I +C I D 1. R ' D I. C 1 

-IT T 

-(A-BR dJ.Cj) 

T -1 

c i D i. R 

R X c i 

-1 T 

R B 

R 1 


and therefore the Hamiltonian H^Cy) of (2.2- la) is the A-matrix of r' ! (s). The 
realization has no uncontrollable or unobservable modes on the jco-axis. Thus, HL(y) 
has no eigenvalues on the jco-axis if and only if T' 1 (s) has no poles on the jco-axis. 
Hence, a x , the infimum of y such that ELfy) has no jco-axis eigenvalues, is equivalent 
to the supremum of y such that 


det [T(ja>)] = 0 for all <o. 
That is, a x is the supremum of y such that 


G^Oco^nOco) G* n 0'co)G 12 (jco) 


V 


Ant 

0 

V 

G* 2 (jco)G u 0'a>) G* 2 (jco)G 12 (jco) 


A 


_ 0 

0 _ 

A 


(2.3-6) 


for all co (2.3-7a) 


or 

G 11 [ 1 - G 12 (G 12 G 12 > 1<3 12 ] G 11 = Y 2 for all co (2.3-7b) 

Therefore, a x is the square root of ^max(G u [I - G 12 (G 12 G J2 ) G 12 ]G n (jco)} 
which completes the proof. 



T -IT 

Remark 2.3-2: A realization of I - G 12 (s)[G 12 (-s)G 12 (s)] G 12 (-s) is given by 


' A - B 2 D nC, 

-b 2 bJ 

B 2 D Ia 

•cMAc, 

-(A-B 2 d[ 2 C,) T 

C>D,2 D L> 

_ -<I-Di 2 Dy 2 )Ci 

D .2 B I 

T 1 

o 

t— ‘ 

K> 

to 


(2.3-8) 


which was obtained from a 4n-th order realization and a mathematical pole-zero 
cancellation. The mathematical pole-zero cancellation was carried out by applying 
similarity transformation and deleting the uncontrollable and/or uncontrollable states 
[ 11 ]. 


Computation of a x , the infimum of y such that HJy) has no jco-axis eigenvalues 
Let 


Tl 2 (co) = ^{Ojjtl-G^CG^G^G^lGjjam)} (2.3-9) 

Then, a x , the infimum of y such that H^fy) has no jco-axis eigenvalues, is the supremum 
of r|(co). There are several efficient algorithms available for searching for the supremum of 
rj(co) [8,9]. In [8], a frequency, say co^ is chosen based on the pole location of G n (s) and 
T -IT 

1 * G 12 (s)[G 12 (-s)G 12 (s)] L et Y = Tl(fi>i) and then find all the positive real co's 

such that T) (co) = y. These co's can be easily obtained from computing the jco-axis 
eigenvalues of the Hamiltonian matrix H^y). Now, we have the frequency intervals in 
which Ti(co) > y. Evaluate Ti(co) for each midpoint of these frequency intervals and set y to 
be the maximum of these ri(co)'s. Then, find the new frequency intervals in which T)(co) ^ 
y. According to [9], the convergence of this iterative process is quadratic. This process can 
be repeated until only one frequency interval with rj(co) > yis left and the interval length is 
negligible [9]. In [8], this process is terminated when a frequency interval in which rj(co) is 
greater than y and convex is found. A search method called Brent method was used to 
search for the the supremum of ri(co) in the convex frequency interval. 

Let a Y be the smallest y such that the Hamiltonian matrix J«(y) has no jco-axis 
eigenvalues. The computation of a Y is similar to that of a x . Define 

a = max {a x , a Y } (2.3-10) 

Then a is the infimum of y such that the two algebraic Riccari equations associated to H^y) 
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and J^Cy) have solutions. 

Let y = a and denote the solutions of the two algebraic Riccati equations associated 
to H«»,(y) and J m (y) by X and Y respectively. If X and Y happen to be both positive semi- 
definite and p(XY) < y 2 , then a is the optimal H“ norm according to Theorem 2.2-1. 
Recall that the case in which the optimum occurs at a, the infimum of y such that the two 
Riccati equations have solutions, was referred as case (1). 


Computation of fix, the infimum of y such thatX(y)is positive semidefinite 

Let P x (P Y resp.) be the infimum of y such that X (Y resp.) is positive semi- 
definite. In case (2), Y (or X resp.) is zero for all y > a and hence, the optimum occurs at y 
= P x (or y= Py resp.). In the following, a hyperbolic interpolation search algorithm is 
used to compute Px- 


Denote the eigenvalues of X(y) by ^ ; [X(y)], i=l,2,...,n and let 


f[X(Y)] = 


max{Xj[X(Y)]} if all X.j[X(y)] are positive 
min{Xj[X(Y)] ) if some ^[X^)] is negative 


(2.3-11) 


2 2 

We observe that f[X(y)] is a hyperbola-like function of y 2 . When y 2 equals to P x +£ (or P x - 
e, resp.) and e approaches to zero, f[X(Y)] will approach to infinity (or negative infinity, 

2 

resp.). As y 2 increases from P x +£ to infinity, f[X(y)] decreases monotonically to a positive 
number. 


To search for Px, we will use two passes in our algorithm. The first pass is used to 

r— 2 

find a lower bound x L such that exists and xl < Px and an upper bound xy which 

a 2 a 2 

is greater than P x . The second pass is to search for P x . 

Since the shape of the graph ( y 2 , f[X(Y)]) resembles that of a hyperbola, a 
hyperbolic interpolation search algorithm will give fast convergence in computing Px- For 
notational simplicity, y 2 and f[X(y)] are denoted by x and y respectively in the following. 
In the algorithm, we start from arbitrary three points on the graph (x,y). It is easy to choose 
three positive numbers x ; , i=l,2,3, with x 1 <X 2 <X 3 , such that their corresponding Riccati 

equations associated to FLC/x") have positive semidefinite solutions X(/x”). Computing 
y ; = f[X(7T)] from (3-11), we have the trio (x^), i=l,2,3. Obviously, P x is smaller 
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than and hence x t qualifies as an upper bound of P^, denoted by xy. Let (xy.yy) = 
(xj.yi). The abscissa of any point (x L ,y L ) on the graph (x,y) with y L <0 can serve as a 
2 

lower bound for P x . 

Now, a hyperbola which interpolates these three points, (x^y^, i=l,2,3, can be 
easily determined. That is, the three parameters a, b, and c in the hyperbolic equation, 

(x-a) (y-b) = c (2.3-12a) 

are determined by (Xj,y;), i=l,2,3. They are, 


a 



b=[ 


x 2 y 2 

XiYi 


x, 

1 --^- 


0-f)a 

X 9 


yi 


(2.3- 12b) 


(2.3- 12c) 


c = (x 1 -a)(y 1 -b) (2.3-1 2d) 

In the graph of the hyperbola, y ±<» as x — > a. If the trio are in the neighborhood of 

R 2 O 2 

(p x ,±<*>), the a of (2.3-12b) will be closer than any to P^ Let x a =a and evaluate 
y a =f[X(/x^)] if X(/xj exists. If y a >0, we update (xy.yy) and the trio as follows, 

( x u>yu) = (x a .y a )- ( x 3»y3) = (x 2 .y 2 ). ( x 2 ,y 2 ) = (Xi.yi), and (x^) = (x a ,y a ). 

2 

If y a <0, X a qualifies as a lower bound for P x and therefore we have (x L ,y L ) = (x a ,y a ). If 

X(Jxl) does not exist, i.e., H^/x^) has eigenvalues on the jto-axis, we suggest to find a x 
first. a x , the infimum of y such that Ffjy) has no jOD-axis eigenvalues, can be easily 
computed by using Theorem 2.3-1 and the algorithm in [8,9]. If X(a x ) is positive 
semidefmite, then we have case (1) and the optimal H°° norm is a x . Otherwise, we have 

case (2) and the optimum occurs at Px which still remains to be found. Let x L =a£ and then 

2 

x L qualifies as a lower bound for p x since y L =f[X(/\)] i s negative. 
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Once (x L ,y L ) and (xy.yu) are available, we are ready to go to the second pass to 
2 

search for P x . During the first pass, it is possible to have y a >0 all the time and therefore 
there is no chance to obtain x L . In this event, we do not need the second pass. Instead, we 
just keep updating the trio and the point (x a ,y a ) is moving closer and closer to the optimal 
2 

point (P x ,±oo). The iterative process can be terminated when ly a l is large enough. 

In the second pass, we start with the new trio (x L ,y L ), (xy.yu) and (x E ,y E ). In 
addition to (x L ,y L ) and (xy.yy), the point (x E ,y E ) is the best point which can be chosen 

from the trio at the end of the first pass. The best point here means the one with the largest 
lyjl. Let (Xj.yj) = ( x L *y L ), ( x 2 ,y 2 ) = 0^^). and (x 3 ,y 3 ) = ( x E ,y E ), then a hyperbola 

which interpolates these three points can be obtained from the formulas in (2.3-12a) - (2.3- 

12d). If x L < a < Xy, let x G =a. Otherwise, let x G -( x L +x (J )/2. Evaluate y G =f[X(/x^)]. The 
point ( x L ,y L ) (or (x^y^, resp.) will be updated by ( x G ,y G ) if y G <0 (or y G >0, resp.). 
Before the updating of ( x L ,y L ) (or (x^y^, resp.), we like to update (x E ,y E ) by ( x L ,y L ) 
(( x u>yu) res P-) if lye 1 < tyj (or IyE* < lyuU resp.). This iterative process is repeated until ly G l 
is sufficiently large or the gap Ixy-xj is small enough. In the second pass, we use hybrid 

search scheme: the hyperbolic interpolation to speed up the convergence and the bisection 
to guarantee the convergence when the hyperbolic interpolation does not work. Once the 
trio reaches the neighborhood of the optimum, the convergence rate is quadratic. 

Computation of the lower and upper bounds for the optimum of Case (3 ) 

As mentioned, case (3) is the most likely one to happen. In case (3), the optimum 
occurs at the y such that p[X(y)Y(y)] = y2 where p(XY) is the spectral radius of XY. In 

subsection 2.4, a fast iterative algorithm is proposed to compute the optimum. To be able to 
use that algorithm, initial lower and upper bounds, y L and y v , are needed. Any y with 

which both X(y) and Y(y) are positive semi-definite and p[X(y)Y(y)] > y 2 qualifies as a 
lower bound. Similarly, any y with which both X(y) and Y(y) are positive semi-definite and 
p[X(y)Y(y)] < y 2 can serve as an upper bound. It is easy to have an upper bound which can 
be chosen arbitrarily large. However, choosing a lower bound is not trivial at all. A small y 
may not satisfy the positive semi-definiteness condition and a larger y may violate the 
condition p[X(y)Y(y)] > y 2 . 

Recall that P x (Py resp.) is the infimum of y such that X (Y resp.) is positive semi- 
definite. Let 
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P = max (P x , P Y } (2.3-13) 

Then any y > P which satisfies p[X(y)Y(y)] > y 2 qualifies as a lower bound. In order to 
find a lower bound, it is not necessary to compute Px and P Y although they can be obtained 
from the hyperbolic interpolation search algorithm described above. An easy way to find a 
lower bound is briefly described in the following. 

Denote the eigenvalues of X(y)Y(y) by X i [X(y)Y(y)], i=l,2,...,n and let 


g[X(y)Y(y)] = 


max { A,j[X(y)Y(y)] } if all A^XCy) Y(y)] are positive 
min { A.i[X(y) Y (y)] ) if some ^[Xfy) Y (y)] is negative 


We observe that g[X(y)Y(y)] is also a hyperbola-like function of y 2 . Since the shape of the 
graph (y 2 , g[X(y)Y(y)]) resembles that of a hyperbola, a hyperbolic interpolation search 
algorithm will give fast convergence in finding a lower bound, y L , which satisfies y L > P 
and g[X(y L )Y(y L )] > y L 2 . In the following, y 2 and g[X(y)Y(y)] are replaced by x and y 
respectively to simplify the notation. To find y L , we start from arbitrary three points on the 
graph (x,y). It is easy to choose three positive numbers x i? i= 1,2,3, with xj<x 2 <x 3 , such 

that their corresponding Riccati equations associated to H^C/xj) (J„ (/*"), resp.) have 

positive semidefinite solutions X(/x7) (Y(/x~), resp.). Computing y t = g[X(,/x^)Y(,/x[)] 
from (2.3-14), we have the trio (x^yj), i=l,2,3. 

The objective is to search for the optimum point (x 0 ,y 0 ) on the graph (x,y) such that 
x 0 =y 0 . In the proposed algorithm, two passes are used to achieve this objective. The first 
pass is to find a lower and an upper bound points (x^yj and (xjj.yu) such that y L >x L and 
x u > yu>0> and therefore x L <x 0 <xu since the graph (x,y) is monotonically decreasing when 
y>0. The second pass is to search for the optimum x 0 . If x 1 <y 1 , we have a lower bound 
point (x L ,y L )=(x 1 ,y 1 ). It is trivial to have an upper bound point (xy.yu) from (x 2 ,y 2 ), or 
(x 3 ,y 3 ), or some other point on the graph with large abscissa. Most of time we have X]>yj, 
and (xu,yu) = (xj.yj). To find a lower bound point (x^.y^), we employ hyperbolic 
interpolations as follows. 

Now, a hyperbola which interpolates these three points, (xj.yj), i=l,2,3, can be 
easily determined. That is, the three parameters a, b, and c in the hyperbolic equation of 
(2.3-12a) can be evaluated by (2.3-12b), (2.3-12c), and (2.3-12d) respectively. Let x P be 
the abscissa of the intersection point of the hyperbola and the line y=x and let x a =a. 

Evaluate y a =g[X(/^)Y (JT a )] and y ? =g[X(J^)Y(J^,)] if they exist. We may have a 
lower bound point (x L ,y L )=(x P ,y P ) if y P >x P or have (x L ,yL)=(x a ,y a ) if y a >x a . If (x L ,y L ) is 
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still not available, we will update (x^yy) and the trio by the information of (x P ,y P ) and 
(x a ,y a ) and repeat the iterative process until (x L ,y L ) is found. If y P does not exist (and 
neither is y P ), we will firstly compute a, the infimum of y such that H„,(y) and J^fy) have 
no jco-axis eigenvalues. Recall that a = max {a x , ot y } and a x and a Y can be easily 
obtained by using Theorem 2.3-1 and the algorithms in [8,9]. For any x > a 2 , 
y=g[X(/x )Y (/x )] always exists. Denote a lower bound for the lower bound x L by x LL 
which is set to be a at the first time y P does not exist. When the hyperbolic interpolation 
search does not work (i.e. when x P <x LL ), we let x B = (x LL +Xu)/2 and compute 

y B =g[X(y^)Y(Vv]. If y B >x B , we have (x L ,y L )=(x B ,y B ). Otherwise, either x LL or x v 
will be updated by x B and the trio can be also updated by using (x B ,y B ). 

When a lower and an upper bound points (x^l) and (x^yy) are available, the 
second pass of the algorithm can be used to search for the optimum x 0 . The second pass of 
the algorithm will be stated in the following section. 

2.4 Computation of the Optimal H~ Norm 

Recall that the optimal H°° Norm is the smallest ysuch that the two Riccati equations 
associated to H^Xy) and J^Cy) have positive semi-definite solutions X(y) and Y(y) and 
p[X(y)Y(y)] < y 2 . There are three cases to consider. The optimum may occur at y = a 
which is referred as case (1) and may also occur at either y = Px or y = P Y which is referred 
as case (2). a, P x , and P Y were defined in the previous section. In case (3), the optimal H~ 

Norm is the y o such that p[X(y o )Y(y o )] equals to Y^. 

To compute the optimum, we start from any three points on the graph (y 2 , 
g[X(y)Y(y)]) where g[X(y)Y(y)]) was defined in (2.3-14). If Y(y) or X(y) is zero and the 
optimum occurs at either y = p x or y = p Y , we have case (2). The computation of P x (or 
P Y ) has been addressed in the previous section. The optimum can also occur at y = a x or y 
= a Y , which was referred as case (1). a x or a Y can be easily computed by using Theorem 
2.3-1 and the algorithms in [8,9]. For the rest of the section, we assume that Y(y) and X(y) 
are both nonzero and there are only cases (1) and (3) yet to be considered. In case (1), the 
optimum occurs at y = a and a = max{a x ,a Y ] which can be easily obtained as mentioned. 

If we can find a lower bound y L such that g[X(y L )Y(y L )] > then we have case (3). 
Unless such a lower bound y L does not exist, the computation of a is not necessary. 

For notational simplicity, we replace (y 2 , g[X(y)Y(y)]) by (x, y). The starting three 
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points on the graph (x, y) are (x^y;), i= 1 ,2,3, with x 1 <x 2 <x 3 . First of all, we will try to 

find a lower bound x L (i.e., and an upper bound xu (i.e., 7^). Usually, the Xj's of the 
initial trio (x^y;), i=l,2,3, are chosen as large numbers to guarantee solution existence for 
the Riccati equations. An upper bound x y can be initially assigned as Xj if y t is positive 
and y t < x } . The locating of a lower bound was briefly described at the end of the previous 
section. Recall that any point (x L ,y L ) on the graph (x,y) which is to the right of (3 and 
satisfies x L <yL qualifies as a lower bound point. A hyperbola interpolating these three 
points, (Xj.yj), i=l,2,3, can be easily obtained from (2.3-12). Let x P be the abscissa of the 
intersection point of the hyperbola and the line y=x and let x a be the abscissa of the singular 
point of the hyperbola, i.e., x a =a where a is from (2.3- 12b). These two points (x a ,y a ) and 
(x P ,y P ) are good guesses for (x^y^. If y P £ x P , we have (x L ,y L )=(x P ,y p ). Otherwise, try 
to see if we have y a > x a which gives (xL,yL)=(x a ,y a ). If the above two inequalities are not 
satisfied, we will use (Xg.yJ and (x P ,y P ) to update (xy.yy) and the trio and repeat the above 
iterative process until a lower bound point (x^yt) is procured. 

In the above hyperbolic interpolation search process, y P =g[X(^/x^)Y( 1 /x^)] may 

not exist since the two Riccati equations may not have solutions. If either KLC/x^) or 

J.oC/^p) has eigenvalues on the imaginary axis, we may have case (1). The infimum of y 
such that H„(y) and J^fy) have no eigenvalues on the imaginary axis, i.e., a, should be 
computed by the algorithm proposed in the previous section. If the Riccati solutions X(a) 
and Y(a) are both positive semi-definite and p[X(a)Y(a)] < a 2 , then a is the optimal H~ 
norm. If either X(a) or Y(a) is indefinite, a bisection search need to be used together with 
the hyperbolic interpolation to guarantee convergence. The detailed algorithm is listed at the 
end of this subsection. 

Now, the lower and upper bound points, (x L ,y L ) and (xy.yy) are available. It is 
well known that in the graph (x,y), y is a monotonically decreasing function of x when x > 
P 2 . From the observation of numerous examples, we have a conjecture that y is a convex 
function of x when x > p 2 . Furthermore, the shape of the graph (x,y) resembles a 
hyperbola and therefore a hyperbolic interpolation search method can be used to speed up 
the convergence. Although the proposed algorithm is devised by exploiting the conjectural 
properties of hyperbola resemblance and convexity, it still works even these properties fail 
to hold. In case that any of these properties does not exist, the algorithm will detect it and 
switch the search scheme to the bisection and the convergence will still be guaranteed. Even 
the conjecture is not absolutely correct, the chance for it to fail in practice would be very 



rare. 


To search for the optimal H°° norm, y , which square is the abscissa of the point on 
the graph (x,y) with x=y, we start from a new trio (x L ,yL), (xu,yu), and (x E ,y E ) which are 
available at the end of the first pass, i.e., the pass of finding lower and upper bounds x L 
and xjj. The point (x E ,y E ) is the available one which has the largest ly^ excluding (x L ,y L ) 
and (xy.yu)- We also initialize &l =x l an< i x u =x u- A hyperbola which interpolates these 
three points can be easily determined. Let (x^) = (x L ,yL>, (x 2 ,y 2 ) = (xu>yu)< and (x 3 ,y 3 ) 
= (x E ,y E ), then a hyperbola which interpolates these three points can be obtained from the 
formulas in (2.3- 12a) - (2.3- 12d). The abscissa of intersection point of the hyperbola and 

the straight line y=x, denoted by x G , will be close to x q = 7^. According to our conjecture, 
x G is supposed to be between & L and fcy. If that is not the case, x Q is reassigned as 

(fcj+fcy)#. Now, we evaluate y G = p[x(^x^ )Y(^x^ )]. 

In the graph (x,y), if y is a convex function of x, the following theorem can be 
employed to find a bracket which encloses the optimum x q . 

Theorem 2.4-1: 

Given three points ( x L ,y L )> (Xy.yy), and (x Q ,y G ) on the graph (x,y). Assume that y L > 
x L , y G < x G , and X L < x Q < x G . If y is a convex function of x, we have the following: 


(i) If y G < x Q , then 

06 < x 0 < e 5 

(2.4- la) 

(ii) If y G >x G , then 

05 < x o < 06 

(2.4- lb) 

where 

0 5 

y E ■ 

_ yh ' Xlx l' x g 

(2.4-2a) 

1 

1 

and 

0 6 

X L’ X G 

y -x yu ' y ° 
_ y ° Gx u- x g 

j yu-y G 

(2.4-2b) 


c 

1 

X 

0 


Proof: It is straightforward to show these inequalities by geometry. 
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(2.4-3a) 


Now we have a bracket which encloses the optimum x q , i.e., 

& r < x < X., 

L o u 

where 

£ L = 0 6 and ^ = 85 if y G < X G > (2.4-3b) 

^ l = 0 5 and ^ = 06 if y G > x G - (2.4-3c) 

If the gap between $ L and ^ is negligible, the search algorithm will be terminated and then 
we have the optimum x q = Xy. Otherwise, the trio ( x L ,y L ), (Xy.yy), and (x E ,y E ) shall be 
updated and the hyperbolic interpolation search shall be repeated until the gap between x L 
and is small enough. The updating of the trio is executed as follows. When y G > x G , x G 
is a better lower bound than x L and therefore we update (x E ,y E ) by (x L ,y L ) if lx G -x L l < lx Q - 
x E l and update (x L ,y L ) by (x G ,y G ). When y G < x G , x Q is a better upper bound than x G and 
therefore we update (x E ,y E ) by (xy.y^ if Ixq-xJ < IxQ-Xgl and update (xy.yy) by (x G ,y G ). 

In the above proposed search algorithm, we used the conjectured property of 
convexity. As mentioned above, there is no proof yet for this conjecture. However, it is 
easy to check if the final approximately equals to x q . Evaluate ^ = 

p[X(^ )YC check if ^ ^ and (& G - P v ) is negligible. If the condition is 

satisfied, practically we have the optimum x q =^ u , i.e., we have the optimal H“ norm 

^o = V^u' ^ * s impossible for the condition to fail. Nevertheless, we can update the 

lower and upper bounds x L and x G and use the bisection search to reduce x G -x L in case that 
the condition fails. 

To summarize, an iterative algorithm for computing the optimal H~ norm is listed as 
follows. 

Iterative Algorithm for Computing the Optimal H~ Norm 

Initialization: 

A 1 . Initialize the tolerances £j, e 2 , e 3 , which are small positive numbers, and let x LL = 
0 , x L = 0 and Xy = 0 . 

A2. Arbitrarily choose three large x., i=l,2,3, such that x 3 > x 2 > Xj > 0 and their 
corresponding Riccati solutions XC/x7 ) and Y(/x[ ) are all positive semi-definite. 
A3. If both X and Y are nonzero, go to Dl. 
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Computation of fa, the optimum for Case (2): (First Pass) 

B 1 . Either X or Y is zero. Without loss of generality, we assume that Y=0 and let y. = 
royxT)], i=l,2,3, where f[X(/x[)] was defined by (2.3-11). Let (x^yy) = 

(Xj.yj). 

B 2 . Let x = a where a is evaluated from (2.3-12b). Compute y = f[X(/xJ] if X(/x^) 
exists. If X(/xJ does not exist, go to Cl. 

B3. If y a > (l/zf), the optimal H“ norm is J3Q. Stop. 

B 4 . If y a > 0, update (Xy.yy) and the trio as follows: 

( x u,yu)=(x a ,y a ), (x 3 ,y 3 )=(x 2 ,y 2 ), (x 2 ,y 2 )=(x 1 ,y 1 ), (x 1 ,y 1 )=(x a ,y a ), Go to B2. 

B 5 . (x L ,y L )=(x a ,y a ), (x E ,y E )=(x 2 ,y 2 ). 

Computation of fa, the optimum for Case (2): (Second Pass) 

B6. Let (x 1 ,y 1 )=(x L ,y L ), (x 2 ,y 2 )=(x u ,y u ), (x 3 ,y 3 )=(x E ,y E ), and evaluate the a of (2.3- 
12b). 

B 7 . If x L < a < Xy, let x G = a and go to B9. 

B8. Let x G = (Xj^+Xy)^. 

B 9 . Evaluate y Q = f[X(-/*G)] . 

BIO. If y G > (l/Gj), the optimal H“ norm is Stop. 

B 1 1 . If Ixy - xj < e , the optimal H~ norm is Stop. 

B 1 2 . If y G < 0, then update (x E ,y E ) by ( x L ,y L ) only if ly L l > ly E l, 
and update ( x L ,y L ) by (x G ,y G ) and go to B6. 

B 1 3 . Update (x E ,y E ) by (x u ,y u ) only if ly^jl > ly E l, 
and update (Xy.yy) by (x G ,y Q ) and go to B6. 

Computation of a^, the optimum for Case (1) with Y=0: 

Cl. Compute c^, the infimum of y such that H^fy) has no joo-axis eigenvalues, by 

using Theorem 2.3-1 and the algorithm in [8,9]. 

C 2 . If X(<x x ) is positive semidefinite, the optimal H~ norm is c^. Stop. 

C 3 . Let x L = a£ and and evaluate y L = f[X(J\)] . 

Let (x £ ,y E ) = (x 2 ,y 2 ) and go to B6. 

Find upper and lower bounds for Case (3 ): 

D 1 . Compute y. = g[X(^)Y(V*I>], i=l,2,3, where g[X(y)Y(y)] was defined by (2.3- 
14). 
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D2. If x^yj, then let (* u ,y u )=(x 1 ,y 1 ) and go to D7. 

D3. Let (x L ,y L ) = (x^yp. 

D4 . If X 2 >y 2 , then let (x u ,y u )=(x 2 ,y 2 ) and go to FI. 

D5. If x 3 >y 3 , then let (x u ,y u )=(x 3 ,y 3 ) and go to FI. 

D6 . Double and evaluate y 3 = g[X(/xJ)Y C/^J)] and go to D5. 

D7. Evaluate a, b, c from (2.3- 12b) - (2.3-12d). 

Let x a = a and x p = y [ a+b+V (a+b) 2 - 4(ab-c) ]. 

If x p < x LL , go to D18. 

D 8 . If (x LL =0) and (either X(./x^) or Y (-/*p) does not exist), then go to El. 

D9. Evaluate y p = g[X(/^)Y(/V)]. 

If y P > x p , then (x L ,y L ) = (x p ,y p ), and go to FI. 

DIO. If ( X LL =0) and (x p > y p > 0), then update the following: 

(Xu,yu)=(Xp,y p ), (x 3 ,y 3 )=(x 2 ,y 2 ), (x 2 ,y 2 )=(x 1 ,y 1 ), (x 1 ,y 1 )=(x p ,y p ), 
and go to D13. 

Dll. If (x LL *0) and (x p > y p > 0), then update the following: 

(x u ,y u )=(Xp,yp), (x 3 ,y 3 )=(x 2 ,y 2 ), (x 2 ,y 2 )=(x p ,y p ), 
and go to D13. 

D 1 2 * ( x LL’yLi>( x P’y P )> ( x i-yi)=(x p ,y p ), and go to D7. 

D 1 3 . If (x ll <x ll ) or (x a does not exist), go to D7. 

D 1 4 . Evaluate y a = g[X(/xI)Y (/xj] . 

If y a > X a* 111611 ow = ( x a ’y a ^ and §° to F1 - 

D 1 5 . If (x LL =0) and (x a > y a > 0), then update the following: 

( x U’yu)-< x a >y a )’ (x 3 ,y 3 )-(x 2 ,y 2 ), (x 2 ,y 2 )— (Xj,yj), (Xj,yj)=(x a ,y a ), 
and go to D7. 

D 16. If (Xjj^O) and (x a > y a > 0), then update the following: 

( x u ,y u )=( x a >y a ). (x 3 ,y 3 )=(x 2 ,y 2 ), (x 2 ,y 2 )=(x a ,y a ), 
and go to D7. 

D17. (x LL ,y LL )=(x a ,y a ), (x 1 ,y 1 )=(x a ,y a ), and go to D7. 

D 1 8 . Let x B = (xll+x^/ 2 and evaluate y B = g[X(/xJj)Y (-/*^)] . 

If yB > x b’ then ( x L ’y L ) = ( x B’yB^ md s° to F1 

D 19. If x B > y B > 0, then update the following: 

( x u’yu) = ( x B’y B )’ ( x 3 >y 3 )=(x 2 ,y 2 ), (x 2 *y 2 ) = (x B ,y B ), 
and go to D7. 
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D20. (x LL ,y LL )=(x B ,y B ), (x 1 ,y 1 )=(x B ,y B ), and go to D7. 

Computation of a, the optimum for Case ( 1 ): 

El. Compute a, the infimum of y such that H^y) and J^fy) have no jco-axis 
eigenvalues, by using Theorem 2.3-1 and the algorithm in [8,9], 

E2. If X(a) and Y(a) are positive semidefinite, and g[X(a)Y(a)] < a 2 , the optimal H~ 
norm is a. Stop. 

E3. Let ( x LL ,y LL ) = (a 2 ,g[X(a)Y(a)]) and go to D18. 

Search for the optimum x D at which y 0 -x 0 : 

F 1 . Now we have the trio, ( x L ,y L ), (Xy.yy), and (x E ,y E ), where (x E ,y E ) is one of the 
(x.,y.)'s other than (x L ,y L ) and (Xy,yy) which has the largest ly.l. 

F 2 . Let & L = x L , and x y = Xy. 

F3 . If I Xy - I < Bj, the optimum x 0 = £y . Stop. 

F4. Let (x 1 ,y 1 )=(x L ,y L ), (x 2 ,y 2 )=(Xy,yy), and (x 3 ,y 3 )=(x E ,y E ). Compute a, b, and c 
from (2. 3- 12b) - (2.3-12d). 

F5. x G = y [ a+b+7 (a+b) 2 - 4(ab-c) ], yo = XJXf^Y 
F 6 . If & L < x Q ^ &y , Go to F8. 

F 7 . Xq = (x L +$y)/2, y G = Xj [X(Jx^)Y(Jxq)] . 

F8. If y G < x G , Go to F10. 

F9. a L = max{0 5 ,£ L }, £y = min{0 6 , £y}. Go to Fll. 

F10. = max{0 fi , ^ L }, ^y = min{0 5 , ^y). 

Fll. If IXy-^eygotoFtt. 

FI 2. If 0 < &y - p[X(JV v )Y(jr v )] < e 2 , then the optimum x 0 = fty. Stop. 

F 1 3. Let x G = (x L +Xy)/2, y G = 

FI 4. If y G > x G , then let (x E ,y E ) = ( x L ,y L ), ( x L ,y L ) = ( x G ,y G ) and go to F2. 

F15. Let (x E ,y E ) — (Xy,yy), (Xy,yy) = (x G ,y G ) and go to F2. 

F 1 6. If y G > x G , then let (x £ ,y E ) = (x L ,y L ), ( x L ,y L ) = (x G ,y G ) and go to F4. 

F 1 7 . Let (x E ,y £ ) = (Xy,yy), (Xy.yy) = (x Q ,y G ) and go to F4. 
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2.5. Construction of Optimal Controllers 

From the Glover and Doyle's formulas in Theorem 2.2-1, a suboptimal H°° 
controller can be easily constructed. However, as y approaches to the optimum we will 
encounter the inversion of a singular matrix except case (1) which seldom occurs. To 
eliminate the numerical difficulty, Glover and Limebeer et. al. [10] rederived the optimal 
controller formulas in a descriptor form (or generalized state-space representation). 


Construction of Optimal H°° Controllers for Case (3): 

The Glover and Doyle's formulas in (2.2-3a) - (2.2-3h) can also be written in a 
descriptor form after slight rearrangement. First of all, we consider case (3) which occurs 
much more often than the other two cases. When y reaches the optimum, y 0 , which 

satisfies To = p[X(y 0 )Y(y 0 )], the matrix Z in (2.2-3f) will become infinity since the matrix 

I-^Y(y 0 )X(y 0 ) is singular. If we try to apply the formulas (2.2-3a) - (2.2-3h) directly to 
construct an optimal H°° controller, a numerical difficulty will arise in the implementation of 
the A and fc matrices. We will rearrange these formulas such that an optimal H“ controller 
can be constructed without any numerical difficulty. 

The dual system of the realization in (2.2-3a) can be easily rewritten in a descriptor 
form. The state equation (generalized state equation) of the descriptor representation can be 
split into two set of equations: one involves first derivative of some state variables and the 
other is just an algebraic equation. The state variables which have no derivative in the 
equations can be eliminated and then we have a lower order state space representation for 
the dual system. The dual of the dual system is identical to the original and therefore we 
have an optimal H" controller as follows: 


' A c 

B c ‘ 

. C c 

°c . 


where 

a c = [ v^AdU! - vXu^Adu/v^AdU, ] 

B c = V]B D - V[A D U 2 (VjA D U 2 ) t V% 

C c = [ c D U! - c D u 2 (V^A D U 2 ) + v t 2 a d jj 1 ] z; 1 


(2.5-1) 


(2.5-2a) 

(2.5-2b) 

(2.5-2c) 
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and 


D c = 6 -c d u 2 (vJa d u 2 ) + vX 


(2.5-2d) 


B d = -H 2 + (B 2 + H 12 )ft 

Cq = F 2 " ft(C 2 + F 12 ) 

A d = (B 2 + H 12 )C d + (A + HC)E d 


„-2 


E d = I-7oX(y 0 )Y(y 0 ) 

Lj, Uj, U 2 , Vj, and V 2 are obtained from the 


u d 


'i 


0 


n 


Lo oJ 


(2.5-3a) 

(2.5-3b) 

(2.5-3c) 

(2.5-3d) 

singular value decomposition of E D , i.e., 

(2.5-4) 


Construction of Optimal H°° Controllers for Case (2): 

In case (2), either Y or X is zero and then the optimum occurs at the smallest y such 
that X(y) or Y(y) is positive semidefinite. In the following, with loss of generality we 
assume Y = 0. The optimum occurs at Y = Px> i- e -> the smallest y such that X(y) is positive 
semidefinite. If we try to apply the formulas (2.2-3a) - (2.2-3h) directly to construct an 
optimal H°° controller, a numerical difficulty will arise since X(Px) is infinity and so are the 
matrices A and £. 

Let X.lHooCPx)] he an n-dimensional spectral space corresponding to eigenvalues in 
Re s < 0. Partitioning the matrix constructed from the basis vectors of %_[EL(P X )], we 
have 


%.[IUp x )] = Im 



(2.5-5) 


where Xj, X 2 € IR . Xj is singular and therefore X = X 2 X^ is infinity. To eliminate the 
numerical difficulty, we do the singular value decomposition for x{ as follows, 




U 2 ] 


Lo 



rn 


V 1 


_ _T 


V 

v 2 


(2.5-6) 


25 



Employing the same techniques, the duality and the elimination of unnecessary state 
variables, which we just used for case (3), an optimal H~ controller can be constructed as 
follows without any numerical difficulty. 


VC*) 



(2.5-7) 


where 


A<; = [ vIAdUj - v]a d U 2 (V^A d U 2 ) + V^AdU, ] Z'/ 

(2.5-8a) 

b c = v]b d - v]a d u 2 (v^AdU^ vj b d 

(2.5-8b) 

c c = [ Colli - C D U 2 (V 2 A d U 2 ) + V 2 A D Ui ] z; 1 

(2.5-8c) 

D c = 5 - C D U 2 (V^AoUz) 1 VjBo 

(2.5-8d) 

b d = -h 2 + (b 2 + h 12 )6 

(2.5-9a) 

C D = P 2 - 5(C 2 x 1+ f 12 ) 

(2.5-9b) 

A d = (B 2 + H 12 ) C D + (A + HC) Xi 

(2.5-9c) 


Construction of Optimal H°° Controllers for Case (1): 

For case (1), the formulas (2.2-3a) - (2.2-3h) can be used to construct an optimal 
H~ controller without any numerical difficulty. 


2.6. Illustrative Examples 

Four illustrative examples will be included in this section. Example 1 is a simple 
four-block H°° optimization problem which belongs to case (3). Example 2 is a two-block 
H~ optimization problem which also belongs to case (3). Example 3 is a four-block H°° 
problem which belongs to case (1). Example 4 is a two-block H“ problem which is used to 
illustrate case (2). 
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Example 1: 


The following is a simple four-block H°° optimization problem which is used to 
illustrate the proposed iterative algorithm of computing the optimal H°° norm. A realization 
of the plant G(s) is given by 
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Two passes are used to compute the optimum. The first pass is to find a lower and an 
upper bounds which bracket the optimum and the second pass is to search for the optimum. 
We arbitrarily choose three large numbers x^lOO, x 2 =1000, and x 3 = 10000, and evaluate 

= g[X(V^)Y(V^)], i=l,2,3, by (3-14). That is, we start with the trio 

(Xj.yj) = (100,2. 177485978 109461e+01) 

(x 2 ,y 2 ) = (1000,2. 161627203254472e-K)l) 

(X 3 ,y 3 ) = (10000,2.1 6005678 1642087e+01) 

Since x 1 >y 1 >0, qualifies as an upper bound Xy and thus 
( X U’V = (100,2.177485978109461e+01) 

The hyperbola which interpolates the trio can be easily determined by (3-12). The abscissa 
of the intersection point of the hyperbola with the line y=x is 

x g =2.24 1 226047 32098 6e+0 1 . From (3-14), we evaluate y G = and then 

we have, 

(x G ,y G ) = (2.241226047320986e +01, 2.241227 598351088e+01). 

After just one iteration of the first pass, we have 

(x L ,y L ) = (x G ,y G ) = (2.241226047320986e+01,2.241227598351088e+01) 

since x G <y G> Now, we have a new trio which include a lower bound point ( x L ,y L ), an 
upper bound point (Xy.yy), and an additional point (x E ,y E ) which can be chosen from the 
previous trio as, 

(x E ,y E ) = (1000,2. 161627203254472e+01). 
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After just one iteration of the second pass, we find that the optimum x q satisfies the 
following inequality, 

= 2.24122754162167e+01 £ x £ 2.24 1227541 62 17 le+01 = 

L. O U 

Thus, the optimum x q is approximately equal to with accuracy up to 14 digits. For 
double checking, we evaluate 

P[X(V^)Y(V^)] = 2.24 1 227 54 1 62 1 705e+0 1 
which is a little bit less than Therefore, the optimum x 0 = 2.24122754162171e+01. 
That is, the optimal H~ norm is y o = -/*o = 4.734160476390413. 

With y = 4.7341604768 which is very close to the optimum, from (2.2-3) and partial 

fraction expansion we have a suboptimal controller as follows, 

K(s) _ -8.9043458823e+00 . l,0388615552e+ll 

s+8 .7 542343 1 40e-0 1 s+2.1943944664e+10 

The second term of the above expression is a wide band low-pass filter which can be 
approximated by a direct feed-through term when y approaches to the optimum. With 
Y=Y 0 = 4.734160476390413, from (2.5-7) we have an optimal controller as follows 


Ko pt (s) = 


-0.87541981354051831 

-0.1382887384996682 

4.451128374445501 

-4.734160476390407 


which is first-order, one order less than that of the plant 


Example 2: 

The following is a two- block H°° optimization problem which was studied in [12]. 
A realization of the plant G(s) is given by 



In the first pass, we will find a lower and an upper bounds which bracket the optimum. We 
arbitrarily choose three large numbers Xj=10, x 2 =100, and x 3 =1000, and evaluate y. = 

i=l,2,3, by (2.3-14). That is, we start with the trio 
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(x^) = (10,9.570601437485001e-04) 

(x 2 ,y 2 ) = (100,8. 398403602350264e-04) 

(x 3 ,y 3 ) = (1000,8.29555697 1685659e-04) 

Since x^yj^, Xj qualifies as an upper bound Xy and thus 
(x^y^ = (10,9.570601437485001e-04) 

The hyperbola which interpolates the trio can be easily determined by (2.3-12). The 
abscissa of the intersection point of the hyperbola with the line y = x is x Q = 

1 ,237935857327897e+00 and from (2.3-12b) we have x a = a = 
1.237024692928173e+00. From (2.3-14), y G and y a can be computed and then we have 

(x G ,y G ) = (l-237935857327897e+00,4.253228336726583e-02) 

(x a ,y a ) = ( 1 .237024692928 173e+00, 4.39398 1296506005e-02) 

Now, (Xy.yy) and the trio are updated as follows: 

(x u>yu ) = (x a ,y a ) = (1.237024692928 173e+00, 4.39398 1296506005e-02) 

(x 3 ,y 3 ) = (x 1>yi ) = (10,9.570601437485001e-04) 

(x 2 ,y 2 ) = (x G ,y Q ) = (1.237935857327897e+00,4.253228336726583e-02) 

(x r yj) = (x a ,y a ) = (1. 237024692928 173e+00, 4.39398 1296506005e-02) 

After the second iteration of the first pass, we have 

(x G ,y G ) = (1.21 09 871637 64649e+00, 1 . 1 8 22 1 4405325450e+00) 

(x a ,y a ) = (1.21 0025237996426e+00,4. 621151 922644384e+0 1 ) 

Since y a >x a and x G >y G , we can update (x u ,y (J ) and have ( x L ,y L ) as follows: 

(Xjj.yy) = (x G ,y G ) = ( 1 .2 1 0987 1 63764649e+00, 1 . 1 822 1 4405325450e+00) 

(x L ,y L ) = (x a ,y a ) = (1.21 002523799 6426e+00,4.62 1 1 519226443 84e+0 1 ) 

Now, we have a new trio (Xy.yy), ( x L ,y L ), and (x £ ,y E ) to start with in the second pass 
where 

(x E ,y E ) = ( 1 .237024692928 1 73e+00, 4.39398 1296506005e-02). 

After just one iteration of the second pass, we find that the optimum x q satisfies the 
following inequality, 

x = 1.210963712471e+00 £ x n <, 1.21 09637 12497e+00 = 

Thus, the optimum x q is approximately equal to with accuracy up to 12 digits. 
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Example 3: 


The following is a four-block H“ optimization problem which realization is given 
by 
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This realization is almost the same as that of Example 1. The only difference is in the 22 

position of the A matrix. In the first pass, we try to find a lower and an upper bounds 
which bracket the optimum x q . We start with the trio 

(x 1>yi ) = (100,1.879778879610986e-01) 

(x 2 ,y 2 ) = (1000, 1 .873943983436236e-01) 

(x 3 ,y 3 ) = (10000, 1.873362964186965e-01) 

The abscissa of the intersection point of the hyperbola interpolating these three points and 
the line y=x is x G = 5.884887925134553e-01. X(x Q ) does not exist and neither is y Q . By 

Theorem 2.3-1 and the algorithm in [8,9], we compute c^, the infimum of y such that 

H^XY) has no jco-axis eigenvalues, as 

a* = 0.89442719099992. 

It is easy to check that X(a x ) and Y(a x ) are both positive semidefinite and 

0.691300276370028 = [p[X(a x )Y(a x )]] 1/2 < = 0.89442719099992. 

Therefore, = 0.89442719099992 is the optimal H~ norm. 

Example 4: 

The following is a two-block H“ optimization problem which realization is given by 
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The Riccati solution Y is zero for all y. The optimum will either occur at (3 X or a x- Two 
passes are used to search for P x . In the first pass, we will try to use hyperbolic 
interpolations to find a lower bound. If it does not work, a x will be computed. It is easy to 
check if a x is the optimum. If not, a x can be used as a lower bound for P x . Once a lower 
and an upper bounds are available, we go to the second pass to search for P x . In the 
beginning, we arbitrarily choose three large numbers x^lOO, x 2 =110, and x 3 =120, and 

evaluate y. = f[X(7^i)], i=l,2,3, by (2.3-1 1). That is, we start with the trio 

(Xj'.yj) = (100,2.320586064570914e+01) 

(x 2 ,y 2 ) = (110,2.314605879013854e+01) 

(x 3 ,y 3 ) = (120,2.30965287 1 632076e+0 1 ) 

Since y t >0, x t qualifies as an upper bound Xy and thus 

(Xy.yy) = (100,2.320586064570914e+01) 

The hyperbola which interpolates the trio gives 
x = a = 3.560891 3629 13274e+00 

a 

X(/xI) does not exist and neither is y . Now, we compute a x = 5.0000000001. X(a x ) is 
not positive semidefinite, and so a x is not the optimum. However, we can use as a 
lower bound. Let x L =a^ and evaluate y L =f[X(. as follows: 

(x L ,y L ) = (5.000000001, -5.00474676615245 le+00) 

Now, we have a new trio (x^yy), (x L ,y L ), and (x E ,y E ) to start with in the second pass 
where 


(x E ,y £ ) = (x 2 ,y 2 ) = (110,2.3 1460587901 3854e+01). 

After eight iterations of the second pass, we have 

(x a ,y a ) = (2.63 3024605035 1 3 3e+0 1 ,2. 5304960 1621 9245e+0 1 ) 
(x a ,y a ) = (2.509367259661 183e+01, 2.55 148 8425 827425e+01) 
(x a ,y a ) = (2.500069275875799e+01,-2.360866054339716e+01) 
(x a ,y a ) = (2.500 1 1 1983723 194e+01 ,- 1 .35357 1 34805368 1 e+03) 
(x a ,y a ) = (2.5001 1 3497987 190e+0 1 , 1 .906244532943883e+03) 
(x a ,y a ) = (2.5001 12857428468e+01,-1.238558367307747e+05) 
(x a ,y a ) = (2.5001 128671 14679e+01,8.504081686717725e+09) 
(x a ,y a ) = (2.5001 128671 14538e+01,3.624506009602072e+13) 
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We can see that the convergence rate increases rapidly when the trio is in the neighborhood 
of the optimum. The optimal H°° norm is Px = = 5.0001 12865840668. 


3. CONCLUSION AND WORK FOR FUTURE RESEARCH 
3.1 Conclusion 

Glover and Doyle's Riccati equation method is a breakthrough in the solution of H°° 
optimization problems. However, the algorithms available for the computation of the 
optimal H°° norm are slow and inefficient. In this report, a fast iterative algorithm based on 
hyperbolic interpolations and a conjecture of convexity was proposed to compute the 
optimal H°° norm for four-block, two-block, and one-block H°° optimization problems. 
This algorithm is complete which can handle all possible cases and its convergence is 
quadratic when the trio to determine the hyperbola is close to the optimum. The numerical 
difficulty arising in applying Glove and Doyle's formulas to construct an optimal H“ 
controller was also eliminated. 

The H“ tools, the efficient algorithm for computing the optimal H~ norm and the 
reliable construction of optimal H~ controllers, which we have just developed in the 
previous research period will be very useful in our future research. 


3.2 Work for Further Research 

3.2.1 M*A Structure and Robust Stability Analysis for Structured 

Uncertainties 

M-A structure is a rearrangement of a perturbed system where M(s) is the nominal 
system and A is a block diagonal matrix which consists of all the perturbations. M-A 
structure is essential in the SSV analysis and design techniques. Although it is always 
possible to pull out all uncertain parts from a perturbed system and form an M-A structure, 
very little about how to obtain this structure has been addressed in the literature. A 
systematic procedure for constructing M-A structure is proposed. 

A minimal M-A structure means that the dimension of A (or M) of the M-A 
structure is minimal. We can construct an M-A structure for a given perturbed system, but 
the dimension of the structure may be unnecessarily large. A nonminimal structure will 
cause unnecessary complexity in computation and therefore a minimal M-A structure is 
essential in the computation of the SSV. Nevertheless, none about the minimality of M-A 
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structure has been addressed in the literature. A progress in the construction of a minimal 
M-A structure will gready simplify robust stability analysis for structured uncertainties. 


3.2.2 Design of Robust Controllers via H~ Optimization 

The next important issue is how to incorporate robust stability requirement into the 
picture of the optimal controller design. Doyle et. al. [11] included a robust stability 
measure into a cost function together with a robust performance measure and formulated a 
structured singular value minimization problem. The concept of the SSV ball described in 
the previous subsection can be also used to formulate an SSV optimization problem. It is a 
realistic and nonconservative formulation. Nevertheless, no easy solution seems to be 
available for this problem in the near future. Direct incorporating the nonconservative 
robust stability requirement into the controller design may be difficult at present time. 
However, we will keep this problem in mind and consider it as one of our long-term 
research problems. A feasible indirect design approach is proposed. 

In the proposed approach, we will use the H“ norm of the complementary 
sensitivity function as a measure of robust stability and will formulate an H~ optimization 
problem which minimizes the maximum eiror energy subject to a robust stability constraint. 
The error reduction and the robust stability can be traded off by choosing the weighting 
matrices in the cost function. Initially, a weighting matrix is chosen and the corresponding 
H°° optimization problem is solved to obtain an optimal controller. The weighting matrices 
are modified iteratively until the robust stability constraint is just satisfied. The way the 
weighting matrices affect the trade-off is only partially understood. In addition to the 
magnitude of the weighting matrix, the structure of the weighting matrix is also an 
important factor in the trade-off. We will investigate how the weighting matrices affect the 
trade-off. We will also develop an iterative updating procedure for the weighting matrices 
by which a robust controller can be designed such that the closed-loop system is robustly 
stable and the maximum of error energy is minimized. 


3.2.3 Design of Robust Controllers via Mixed H 2 /H~ Optimization 

The optimal H~ controller which minimizes the H~ optimization problem usually 
has wide bandwidth and leads to a poor H 2 performance. In [5], we found that a little bit 
sacrifice of the H°° norm will greatly reduce the bandwidth of the controller and improve the 
H 2 performance tremendously. It is practical to formulate a robust control problem as a that 
of minimizing an H 2 cost function subject to an H~ bound. The H 2 cost function is 
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equivalent to the well known LQG cost function and the H“ bound can take care of robust 
stability and robust error reduction. 

The only approach available in the literature for solving the mixed H 2 /H“ 
optimization problem was proposed by Bernstein and Haddad [13]. A set of coupled 
Riccati equations are involved and only in a special case these coupled Riccati equations can 
be solved. Furthermore, the numerical algorithm for these equations is a homotopy 
algorithm which is not efficient. A research in finding more efficient approach for the 
mixed FP/H" optimization problem is proposed. 


3.2.4 Controllability and Observability of Perturbed Plant 

Although a number of DOC (Degree of Controllability) and DOO (Degree of 
Observability) measures have been defined to quantify the controllability and observability 
of a system, none of them have been used as part of the design of robust controller. 

In this study, we shall address this question. Specifically, first we shall determine 
how sensitivity are the available DOCs and DOOs to the structured uncertainties and 
unmodelled dynamics. Then we shall find a maximal domain of perturbation £ which 
guarantees a specified level of Degree of Controllability (DOC) and another domain of 
perturbation which guarantees a specified level of Degree of Observability (DOO). These 
domains give the controllability and observability margins which take into account the 
uncertainties in the system model. Once these margins have been defined, then they will be 
integrated with the stability margins to develop a methodology for the design of robust 
controllers. 
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